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ABSTRACT 

It has been recently empirically established that some of the directly observed pa- 
rameters of GRBs are correlated with their important intrinsic parameters, like the 
luminosity or the total radiated energy. These correlations were derived, tested and 
used to standardize GRBs, i.e., to derive their luminosity or radiated energy from 
one or more observables, in order to construct an estimated fiducial Hubble diagram, 
assuming that radiation propagates in the standard ACDM cosmological model. We 
extend these analyses by considering more general models of dark energy, and an 
updated data set of high redshift GRBs. We show that the correlation parameters 
only weakly depend on the cosmological model. Moreover we apply a local regression 
technique to estimate, in a model independent way, the distance modulus from the 
recently updated SNIa sample containing 307 SNIa (As tier et al. 2006]) . in order to 
calibrate the GRBs 2D correlations, considering only GRBs with z < 1.4. The derived 
calibration parameters are used to construct a new GRBs Hubble diagram, which we 
call the calibrated GRBs HD. We also compare the estimated and calibrated GRBs 
HDs. It turns out that for the common GRBs they are fully statistically consistent, 
thus indicating that both of them are not affected by any systematic bias induced 
by the different standardizing procedures. We finally apply our methods to calibrate 
95 long GRBs with the well-known Amati relation and construct the estimated and 
calibrated GRBs Hubble diagram that extends to redshifts z ~ 8. Even in this case 
there is consistency between these datasets. This means that the high redshift GRBs 
can be used to test different models of dark energy. We used the calibrated GRBs HD 
to constrain our quintessential cosmological model and derived the likelihood values 
of fl m and u>(0). 

Key words: Gamma Rays : bursts - Cosmology : distance scale - Cosmology : cos- 
mological parameters 



1 INTRODUCTION 

Recent observations of high redshift supernovae of type la (SNIa) revealed that the universe is now expanding at an accelerated 
rate. This surprising result has been independently confirmed by observations of small scale temperature anisotropies of the cos- 
mic microwave background radiation (CMB) (Astier ct al. 2006; Kowalski ct al. 2008; Amanullah ct al. 2010, Jarosik &al 20010). 
It is usually assumed that the observed accelerated expansion is caused by a so called dark energy, with unusual properties. 
The pressure of dark energy pd e is negative and it is related with the positive energy density of dark energy e^e by pd e = wede, 
where the proportionality coefficient w < 0. According to the present day estimates, about 75% of matter-energy in the 
universe is in the form of dark energy, so that now the dark energy is the dominating component in the universe. The nature 
of dark energy is not known. Proposed so far models of dark energy can be divided, at least, into three groups: a) a non zero 
cosmological constant, in this case w = — 1, or b) a potential energy of some, not yet discovered, scalar field, or c) effects 
connected with non homogeneous distribution of matter and averaging procedures. In the last two possibilities, in general, 
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w is not constant and it depends on the redshift z. Observations of type la supernovae and small scale anisotropies of the 
cosmic microwave background radiation are consistent with the assumption that the observed accelerated expansion is due 
to the non zero cosmological constant. However, so far the type la supernovae have been observed only at redshifts z < 2, 
while in order to test if w is changing with redshift it is necessary to use more distant objects. New possibilities opened up 
when it turned out that some of the Gamma Ray Bursts are seen at much higher redshifts, the present record is at z — 8.26 
(Grci ner et al. 2008[) . Neither the type la supernovae nor the GRBs are true "standard candles". A type la supernova is a 
final stage of evolution of a white dwarf, which is a member of a binary system and accreted enough matter to reach the 
Chandrasekhar mass limit of 1.4 Mq. Since the exploding white dwarf always contains mass of about the Chandrasekhar limit, 
and white dwarfs have similar structure, it is reasonable to assume that when they explode they release approximately the 
same amount of energy. Sterling Colgate already in 1979 suggested that type la supernovae could be used as standard candles 
(Colgate 1979). When finally, late in nineteen nineties, astronomers were able to observe very distant type la supernovae, 
and used them to determine distances to high redshift galaxies, they discovered that the expansion of the universe, instead of 
slowing down as expected, is accelerating. Because this conclusion has very important consequences it was prudent to check 
if type la supernovae have the same intrinsic luminosity. It turned out that their intrinsic luminosity varies quite a bit, but 
with the help of several observational characteristics it can be estimated with an accuracy of about 15%. GRBs are even more 
enigmatic objects. First of all the mechanism that is responsible for releasing the incredible amounts of energy that typical 
GRB emits is not yet known (see for instance Meszaros 2006 for a recent review on GRBs). It is also not yet definitely known 
if the energy is emitted isotropically or is beamed. Despite of these difficulties GRBs are promising objects that can be used 
to study the expansion rate of the universe at high redshifts ( |Bradley 2003] ISchaefer 20031 IDai et al. 20041 IBloom et al. 20031 
IFirmani at al. 20051 ISchaefer 20071 ILi et al. 20081 lAmati et al. 20081 ITsutsui et al. 2009J1 . Soon after the BeppoSAX satellite 
was launched in 1996 an afterglow of now famous GRB970228 has been observed and for the first time it was possible to 
determine its redshift and conclusively show that GRBs occur at cosmological distances. Since then, new satellites have been 
launched and it became possible to observationally determine several important parameters of GRBs. First of all, it is now 
possible to determine not only the light curves but also the spectra of GRB radiation for a wide range of photon energies. Using 
the observed spectrum it is possible to derive additional parameters, for example the peak photon energy E pea k, at which the 
burst is the brightest, and the variability parameter V which measures the smoothness of the light curve (for definitions of 
these and other parameters mentioned below, see Schaefer (2007)). From observations of the afterglow it is possible to derive 
another set of parameters, redshift is the most important, and also the jet opening angle Ojet, derived from the achromatic 
break in the light curve of the afterglow, the time lag Ti ag , which measures the time offset between high and low energy GRB 
photons arriving at the detector, and trt - the shortest time over which the light curve increases by half of the peak flux 
of the burst. However the most important parameters - the intrinsic luminosity L, and the total radiated energy _E 7 are not 
directly observable. Assuming that GRBs emit radiation isotropically it is possible to relate the peak luminosity, L, to the 
observed bolometric peak flux Pboio by 

L = 4nd 2 L (z,cp)P bo i , (1) 

where cLl is the luminosity distance, and cp denotes the set of cosmological parameters that specify the background cosmological 
model. Equivalently one can use the total collimation corrected energy _E 7 defined by 

E 1 = 47rd 2 l (z, cp)SboioF beam (l + z) -1 , (2) 

where Sboio is the bolometric fluence and Fbeam = (1 — cos Qjet) is the beaming factor. It is clear from these relations that, in 
order to get the intrinsic luminosity, it is necessary to specify the fiducial cosmological model and its basic parameters. But we 
want to use the observed properties of GRBs to derive the cosmological parameters. The way out of this complicated circular 
situation has been proposed by Schaefer (2007) and others, who studied correlations between the directly observed parameters 
of GRBs (E pe ak, T~i ag , Gjet, and V) and derived parameters like L and Z? 7 . In his detailed study of these correlations Schaefer 
used the standard ACDM and a quintessence model in which w is redshift dependent. He showed that the parameters of 
correlations between L or _E 7 and the observed parameters -E pe afc, Tiag, Qjet, and V only weakly depend on the cosmological 
model. In conclusion his results show that these correlations can be used to standardize the intrinsic luminosity of GRB bursts 
and so to transform them into standard candles. Recently Basilikos and Perivolaropoulos (2008) have extended the work of 
Schaefer, by testing the stability of the correlation parameters with respect to possible evolution with redshift and variations 
of the basic parameters of the background cosmological model. They concluded that the correlation parameters practically do 
not depend on the redshift. They assumed that the background cosmological model is the standard ACDM model but included 
the Q. m parameter into the set of fitting parameters. In fact they obtained very similar values of the correlation parameters 
as those obtained by Schaefer, but the Q m parameter (even at the la level) was poorly determined. Almost at the same time 
and independently Cardone, Cappoziello and Dainotti (2009) extended the analysis performed by Schaefer, by including more 
objects and adding one recently discovered correlation between the GRBs X-ray luminosity Lx and the time T a , which is 
characterizing the late afterglow decay. They then, using a Bayesian procedure, fitted all the correlations, estimating their 
parameters, considering as a fiducial cosmological model the ACDM model, with parameters derived from the WMAP5 data. 
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Finally to check how their results depend on the parameters of the background cosmological model they fitted the correlations 
on a parameter grid assuming five different values for Q, m and allowing the dark energy equation of state parameter w to vary 
with z as w(z) — wo + w a z/(l + z), where (wo, w a ) were restricted to vary in the range — f .3 < too < —0.3 and — f < w a < 1 . 
They have found that the distance modulus defined as 

H{z) = 25 + 5\ogd L {z, cp), (3) 

only very weekly depends on the background cosmological model with variations not exceeding 1%. Using a sophisticated 
statistical procedure and the SNfa data they were able to calibrate the correlation parameters of GRB observables and at 
the same time to determine basic parameters that specify the background cosmological model (Cardonc ct al. 2008). The 
presented so far studies of the correlations between observed parameters of the GRBs used as the background cosmological 
model either the standard ACDM model or a model in which the dark energy equation of state can change by assuming that 
either w(z) is a linear function of the redshift or that w(z) — Wq + w a z/(l + z). Let us note that in recent years a host 
of very different models of dark energy were proposed, with more complicated dependence of w(z) and most of them are 
able to fit at least the SNIa data. It is then very important to investigate how strongly, if at all, the correlation parameters 
depend on the model of dark energy. In particular it is interesting to ask what happens if w(z) is rapidly changing in the 
present epoch and/or when the luminosity distance is noticeably different from the ACDM model only at z > 2. With this 
aim in mind in this note we present results of our analysis of correlations between the observed parameters of the set of 
GRBs assembled by Schaefer and later extended by (ICardone et al. 2008) and derived parameters, first by assuming that 
the background cosmological model is one of the quintessence models and secondly by considering an artificial luminosity 
distance, which at large redshifts gives much larger distances in comparison with the standard ACDM model. To study the 
first problem we selected one of the quintessence models that we have studied some time ago (Rubano & Scudcllaro 2002 
Demianski et al. 2005; Pavlov ct al. 2002; Rubano ct al. 2004), and showed that it is consistent with the basic cosmological 
tests, like type la supernovae, the observed CMB temperature anisotropies power spectrum, and basic parameters of the large 
scale structure as determined by the 2dFGRS and Sloan Digital Sky surveys. In this model the dark energy is represented by 
a self interacting scalar field with exponential potential which is minimally coupled to gravity. The parameters of the potential 
have been chosen in such a way that the coupled Einstein - scalar field dynamical equations describing the flat universe can be 
analytically integrated, what is clearly of a great advantage in our computations. The energy density of dark energy depends 
on time and also the dark energy equation of state parameter w$ is time dependent. In Fig. [T] we show how the u>^> parameter 
is changing with the redshift. As is apparent from this figure w,f, is rapidly changing for < z < 5 and is almost constant and 
equal to —I at earlier epochs. We have picked this model precisely because of this dramatic difference between the behaviour 
of the dark energy in comparison with previously considered background cosmological models. The statistical analysis of the 
correlation relations between the observed parameters of GRBs and their important intrinsic parameters - the total intrinsic 
luminosity L and total emitted energy E 1 performed by Schaefer (2007), Basilikos and Perivolaropoulos (2008) and Cardone 
et al. (2008) shows that the correlation between E pea k and E~/ seems to be the most robust. The total emitted energy is given 
by equation (2) and, as it is apparent from this relation, in order to calculate E-y, it is necessary to specify the background 
cosmological model. We begin our analysis using the quintessense model with exponential potential with parameters fixed by 
fitting this model to the set of type la supernovae data, the power spectrum of CMB temperature anisotropies and parameters 
of the observed large scale structure ( for detailed description see Demianski et al. (2005)). Using the luminosity distance 
derived from this model, we fit our correlations, estimate their parameters, and then we update the estimated GRBs Hubble 
diagram. Then, to check the robustness of the discovered correlations between some of the observed parameters of GRBs and 
their basic intrinsic properties, we have tested the stability of values of the correlation parameters by considering an ad hoc 
definition of the luminosity distance that gives much larger distances to objects at z > 2 than either of the models mentioned 
above. Comparing the likelihood contours in the plane defined by the values of the slope parameter, a, and the intrinsic 
scatter, crjntp, it turns out that our artificial luminosity distance is changing the maximum likelihood values of the correlation 
parameters but the difference is not large. However the situation changes when we consider, in the space of parameters, the 
regions of I a, 2 a and 3 a of confidence for our crazy model. It turns out that with respect to the same regions constructed 
for the quintessential model they overlap only at I a, but differ consistently at higher levels of confidence. 

Moreover we apply a local regression technique to estimate, in a model independent way, the distance modulus from 
the recently updated SNIa sample containing 307 SNIa, in order to calibrate the GRBs 2D correlations considering only 
GRBs with z < 1.4. The derived calibration parameters are used to construct a new calibrated GRBs Hubble diagram. We 
also compare the estimated and calibrated GRBs HDs. We finally apply our methods to calibrate 95 long GRBs with the 
well-known Amati relation and construct the estimated and calibrated GRBs Hubble diagram that extends to redshifts z ~ 8. 
Even in this case there is consistency between these datasets. This means that the high redshift GRBs can be used to test 

1 ""mi is the amount of intrinsic scatter which we expect around the best fit line, due to the empirical nature of the studied correlations, 
as we will clarify in the section [2] 
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Figure 1. The behaviour of the dark energy equation of state parameter on redshift for our fiducial model. It is rapidly changing 
for < z < 5 and is almost constant and equal to —1 at earlier epochs, when the scalar field mimics an effective cosmological constant. 

different models of dark energy. We used the calibrated GRBs HD to constrain our quintessential cosmological model and to 
derive the likelihood values of Q m 
and w(0). 



2 STANDARDIZING THE GRBS AND CONSTRUCTING THE HUBBLE DIAGRAM 



Following (Schacfcr 2007; Cardo ne et al. 20 08 ). we consider six luminosity relations for GRBs. First, we consider five that 
relate GRB peak luminosity, L, or the total burst energy in the gamma rays, E-,, to observables of the light curves and/or 
spectra: ri ag (time lag), V (variability), -E pea k (peak of the vF v spectrum), and trt (minimum rise timefl: 
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where the a; and &; are fitting parameters. To these known correlations we added a new correlation 
log [L x (T a )} = b 6 + a e log [T a /(1 + z)], 



(4) 
(5) 
(6) 
(7) 
(8) 

(9) 



that was recently empirically found by Dainotti et al. (2009) and later confirmed by Ghisellini et al. (2008) and by Yamazaki 
(2008). 

2.1 Fitting the correlations and estimating their parameters 

In fitting these five calibration relations, we need to fit a data array {xt,yi} with uncertainties {(T x ,i,&y,i\, to a straight line 
y = b + ax. (10) 



2 It is worth noting that the -E poa k — E 1 correlation listed below is not physically a 2D correlation, because it involves three obervablcs: 
the measured fluence, the spectral peak energy and the break time of the optical afterglow light curve needed to measure the jet opening 
angle. However it can be treated as a 2D correlation from the statistical point of view (to perform the fit). 



GRBs Hubble diagram 5 



Eg. (|10[) is a linear relation which can be fitted to a given data set (xi, yi) in order to determine the two fit parameters (a, bjj. 
Actually, the situation is not so simple since, both the (y,x) variables are affected by measurement uncertainties (a x ,a y ) 
which can not be neglected. Moreover, a v /y ~ a x /x so that it is impossible to choose as independent variable in the fit 
the one with the smallest relative error. Finally, the correlations we are fitting are mostly empirical and are not yet derived 
from an underlying theoretical model determining the detailed features of the GRBs explosion and afterglow phenomenology. 
Indeed, we do expect a certain amount of intrinsic scatter, aint, around the best fit line which has to be taken into account and 
determine together with (a, b) by the fitting procedure. Different statistical recipes are available to cope with these problems. It 
is not clear how the fitting technique employed could affect the final estimate of the distance modulus for a given GRB from the 
different correlations so that it is highly desirable to fit all of them with the same method. Following Cardone et al. (2008), we 
apply a Bayesian motivated technique (D'Agostini 2005 ) maximizing the likelihood function £(a, 6, a int ) — exp [— L(a, b, a int )] 
with : 

L(a, b, a int ) = - ^ In (af ni + Oy t + a 2 a Xz ) 

_i_ 1 (Vi- ax * ~ *>) 2 , s 

where the sum is over the N objects in the sample. Note that, actually, this maximization is performed in the two parameter 
space (a, ai n t) since b may be estimated analytically (solving the equation — (L(a,b,ai n t)) — 0, as: 



6 = 



^ al 



Vi 



+ a l ai 



V 



+ a^ + aVi, 



(12) 



To quantitatively estimate the goodness of this fit we use the median and root mean square of the best fit residuals, defined 
as 8 = y i,s — yfu which are computed for the different correlations we consider. To quantify the uncertainties of some fit 
parameter p it we evaluate the marginalized likelihood £i(pi) by integrating over the other parameter. The median value for 
the parameter pi is then found by solving : 

£i(pi)dpi = 2 d(pi)dpi . (13) 

tin Jpi,min 

The 68% (95%) confidence range (pij,Pi : h) are then found by solving ([D'Agostini 2005P : 

Ci(pi)dpi = / d(pi)dpi , (14) 

Pi, l J Pi,min 

Ci{pi)dpi = - I £i(p z )dpi , (15) 

i.med Jpi, m in 

with e — 0.68 and e — 0.95 for the 68% and 95% confidence level. Just considering the correlation relation between E y and 
Ep e ak of the form 

log E-, = a log E peak + b , (16) 

we find that the likelihood method gives a — 1.37, b = 50.56 and a^t = 0.25. In Fig. [5] we show the likelihood contours in 
the (a, aint) plane and in Fig. Owe show the correlation between the observed \ogE pea k and derived logi? 7 with our assumed 
background cosmological model. The solid line is the best fit obtained using the D'Agostini's method (D'Agostini 2005) and 
the dashed line is the best fit obtained by the weighted x 2 method. If one marginalizes with respect to b, then the likelihood 
values a and a are a = 1.39 and a — 0.27. In Table 1 we present data for the other correlation relations. It is apparent that 
the correlation parameters only weakly depend on the assumed background cosmological model. Let us stress that in our 
model the dark energy equation of state is changing and the w(z) parameter is quite different from the models considered by 
Schaefer (2007), Basilikos and Perivolaropoulos (2008), and Cardone et al. (2008). To check the robustness of the discovered 
correlations between some of the observed parameters of GRBs and their basic intrinsic properties we have tested the stability 
of values of the correlation parameters by considering an ad hoc definition of luminosity distance that places objects with 
z > 2 at much larger distances than either of the models mentioned above (we call this model a crazy model). First of all in 
Fig. [I] we show the luminosity distance for our quintessence model (solid line) and the crazy model (dashed line). In this case 



3 For the first five correlations we use the sample of GRBs compiled by Schaefer (2007), while the Lx~T a correlation is calibrated 
using a subsample of the Willingale et al. (2007) catalog made out of 28 GRBs with measured redhisft, log Lx(T a ) > 45 and 1 < 
log [T a /(1 + z)] < 5 (see (Dainotti et al. 2008) for the motivation of these limits). Moreover we note that Fi, eam is calculated with the 

— 1/2 

aid of dij{z): for small values of 9j e t, the dependence of F^ eam on di J {z) is found to be: F^ am oc d L . Using the data from Schaefer 
2007 (fi TO = 0.27) we have to multiply the corresponding F^,, am value by the following factor: [d^(0. 27, z)/di(z, CP)] 1 ' 2 
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Figure 2. Regions of 68%, 95% and 99% of confidence in the space of parameters a, <Ti nt . 
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Figure 3. Best fit curves for the E pea k — E^ correlation superimposed on the data. The solid and dashed lines refer to the results 
obtained with the Bayesian and Levemberg - Marquardt estimators respectively. 

using the likelihood minimalization method we get a = 1.69, b = 50.60 and (Tint = 0.30. In Fig. [5] we show the 68%, 95% and 
99% of confidence in the space of parameters a, a n t for our crazy model (red dotted lines), compared with the same regions 
of confidence for the quintessential model: it is apparent that our artificial luminosity distance is changing the values of the 
correlation parameters and the regions of confidence. In Fig. [5] we show the correlation between the observed log E pea k and 
derived log E-, with our assumed luminosity distance function. The statistical analysis of the correlation relations between the 
observed parameters of GRBs and their important intrinsic parameters shows that all the relations have a similar statistical 
weight since both 8 m cd and S rrns have almost the same values over the full set, where 8 — y b s — Ufa- This analysis shows 
only a modest preference for the E-,-E p , L-ri ag and Lx-T a correlations. Taken at face values, the maximum likelihood 
estimates for (Tint favour the E 1 - E p correlation, which seems the most robust, but when one takes into account the 68% 
and 95% confidence ranges all the correlations overlap quite well. In our analysis we have assumed that the fit parameters 
do not change with the redshift, which indeed spans a quite large range (from z — 0.125 up to z = 6.6). The limited number 
of GRBs prevents detailed exploration of the validity of this usually adopted working hypothesis, which we tested somewhat 
investigating if the residuals correlate with the reshift. We have not found any significant correlation. Moreover we tested 
the fit of the £ 7 - E p correlation with respect to the evolution with redshift, separating the GRB samples into four groups 
corresponding to the following redshift bins: z 6 [0,1], z£ [1, 2], z £ [2, 3] and z S [3, 7]. We thus maximized the likelihood in 
each group of redshifts and determined the best fit parameters a, b with la errors and the intrinsic d lSpGrSlOIl (Tint- It turns 
out that no statistical evidence of a dependence of the (a, b, ant) parameters on the redshift exists. This is in agreement with 
what has recently been found by Cardone et al. (2008) and Basilakos & Perivolaropoulos (2008). 

2.2 Constructing the Hubble diagram 

Once the six correlations have been fitted, we can now use them to construct the estimated GRBs Hubble diagram. In order 
to reduce the error and, in a sense, marginalize over the possible systematic biases present in each of the correlations, we take 
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Figure 4. The behaviour of the luminosity distance for our quintessence model (solid line) and the crazy model (dashed line). 
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Figure 5. Regions of 68%, 95% and 99% of confidence in the space of parameters a, (Tint for our crazy model (red dashed lines), compared 
with the same regions of confidence for the quintessential model (solid lines). It turns out that the confidence regions overlap at la, but 
differ consistently at higher levels of confidence. 



a weighted average of the distance modulus provided by each of the six 2D laws in Table 1. Just as an example, let us remind 
that the luminosity distance of a GRB with redshift z may be computed as : 

S 7 (l + z)/4nF beam S bo i , 



d L (z) = < L/47rP 6o ; , 

k L x {l + zf" +2 /47vF x {T a ), 



(17) 




0.0 0.5 

log Epeak 

Figure 6. Best fit curves for the E pf , a ^ — Ej correlation, for our crazy model, superimposed on the data. The solid and dashed lines 
refer to the results obtained with the Bayesian and Levemberg - Marquardt estimators respectively. 
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Figure 7. The calibrated (empty circles) and estimated quintessential (full boxes) Hubble diagram. 



depending on whether some correlation is used. 

The uncertainty of <1l{z) is then estimated through the propagation of the measurement errors on the involved quantities. 
In particular, remembering that all the 2D correlation relations can be written as a linear relation, as in Eq. (|10[) , where y is 
the distance dependent quantity, while x is not, the error on the distance dependent quantity y is estimated as: 



o-(logy) = y/a 2 a 2 (logx) + af nt 



(18) 



and is then added in quadrature to the uncertainties on the other terms entering Eq. (|17[) to get the total uncertainty. The 
distance modulus y,(z) is easily obtained from its definition : 



H{z) = 25 + 51ogd L (z). 



(19) 



with its uncertainty obtained again by error propagation. Following Schaefer (2007) and Cardone et al. (2008), we finally 
estimate the distance modulus for the i - th GRB in the sample at redshift Zi as : 



i 



with the uncertainty given by: 



(20) 



(21) 



where the sum runs over the considered empirical laws. Joining the Willingale et al. (2007) and Schaefer (2007) samples and 
considering that 17 objects are in common, we end up with a catalog of 83 GRBs which we use to build the Hubble diagram 
plotted in Fig. [7] We will refer in the following to this data set as the fiducial GRBs Hubble diagram (hereafter, HD) since 
to compute the distances it relies on the calibration based on the fiducial model. We applied such a procedure first for our 
quintessence model and then for our artificial crazy model. Although the above analysis has shown that the choice of the 
underlying cosmological model has only a modest impact on the final estimate of the distance modulus, we compared our 
estimated fiducial HD with a model independent calibrated HD, carried out using SNIa as distance indicators. Actually the 
SNIa Hubble diagram gives the value of fi(z) for a subset of the GRBs sample with z < 1.4 which can then be used to calibrate 
the correlations (see (|Cardone et al. 2008[l for details). Assuming again that this calibration is redshift independent, one can 
then build up the Hubble diagram at higher redshifts using the calibrated correlations for the remaining GRBs in the sample. 
We implement here a slightly different version of the local regression approach described in Cardone et al. (2008), modifying 
some steps in order to improve the computing time (in particular we use a different choice of the weight function and the 
order of the fitting polynomial-taking ). Once we estimated the distance modulus at redshift z in a model independent way, 
we can now calibrate the GRBs correlations considering only GRBs with z < 1.4 in order to cover the same redshift range 
spanned by the SNIa data. For such subset of GRBs we apply the Bayesian fitting procedure described above to construct a 
new GRBs Hubble diagram that we call the calibrated GRBs HD. Plotted in Fig. [7] this HD now contains only the 69 Schaefer 
GRBs since we have not used the Lx - T a correlation because it is impossible to calibrate it with the local regression based 
method. It turns out that the calibrated and estimated data sets, although containing a different number of objects (83 vs 
69), are consistent (see Figs. [7]). This qualitative agreement is also confirmed by considering the numbers of GRBs that are 
more than 2a away from the theoretical fiducial HD. 
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Figure 8. Regions of 68%, 95% and 99% of confidence in the space of parameters a,ai n t, relatively to the fit of the Amati correlation. 

3 CONSTRUCTING THE HUBBLE DIAGRAM AND FITTING THE E P>I - £iso CORRELATION 

In this section we investigate the possibility of constructing the Hubble diagram from the £" p ,i - -Eiso correlation, here E Pi i 
is the peak photon energy of the intrinsic spectrum and Ei BO the isotropic equivalent radiated energy. This correlation was 
initially discovered on a small sample of BeppoSAX GRBs with known redshifts (Amati et al. 2 002]) and confirmed afterwards 
by Swift observations (Amati 2006). Although it was the first correlation discovered for GRBs it was never used for cosmology 
because of its significant "extrinsic" scatter. However, the recent increase in the efficiency of GRB discoveries combined with 
the fact that E Pt i - Ei so correlation needs only two parameters that are directly inferred from observations (this fact minimizes 
the effects of systematics and increases the number of GRBs that can be used, by a factor ~ 3) makes the use of this correlation 
an interesting tool for cosmology. Previous analyses of the - Ei so plane of GRBs showed that different classes of GRBs 
exhibit different behaviours, and while normal long GRBs and X-Ray Flashes (XRF, i.e. particularly soft bursts) follow the 
E P: i - Ei so correlation, short GRBs and the peculiar very close and sub-energetic GRBs do not (Amati et al. 2008). This 
fact depends on the different emission mechanisms involved in the two classes of events and makes the E P] i - Ei so relation a 
useful tool to distinguish between them (|An tonclli 2009J. A natural explanation for the short /long dichotomy and the different 
locations of these classes of events in the E Pt i - Ei so plane is provided by the "fireshell" model of GRBs (Ruffini 2009). The 
impact of selection and instrumental effects on the E p ,i - Ei so correlation of long GRBs was investigated since 2005, mainly 
based on the large sample of BATSE GRBs with unknown redshifts. Different authors came to different conclusions (see for 
instance (Ghirland a at al. 2005P ). In particular, ( Ghirland a at al. 20 05 ) showed that BATSE events potentially follow the E Ptl 
- E lso correlation and that the question to clarify is if, and how much, its measured dispersion is biased. There were also claims 
that a significant fraction of Swift GRBs is inconsistent with this correlation (IButler et al. 2007]) . However, when considering 
those Swift events with peak energy measured by broad-band instruments like, e.g., Konus-WIND or the Fermi/GUM or 
reported by the BAT team in their catalog ([Sakamoto et al. 2008) it is found that they are all consistent with the E Pt i - 
Ei BO correlation as determined with previous/other instruments (I Amati e t al. 2009]) . In addition, it turns out that the slope 
and normalization of the correlation based on the single data sets provided by GRB detectors with different sensitivities and 
energy bands are very similar. These facts further support the reliability of the correlation ( Amati et al. 2009). In our analysis 
we considered the sample of 95 long GRB/XRF only, compiled in (Amati et al. 2008) and (Amati et al. 2009). The redshift 
distribution covers a broad range of z, from 0.033 to 8.23, thus extending far beyond that of SNIa (z <~1.7), and including 
GRB 092304, the new high-z record holder of Gamma-ray bursts. As above we apply a Bayesian technique (D'Agostini 2005) 
to fit the E Pil - Ei a o correlation, maximizing the likelihood function of Eq. (|11[) . It turns out that our likelihood method gives 
a — 1.56, b — 52.77 and <Ji nt = 0.31. The likelihood contours in the (a, aw) plane are shown in Fig. [8] and finally in Fig. 
Owe show the correlation between the observed \ogE pea k and derived log Ei SO with our assumed background cosmological 
model. The solid line is the best fit obtained using the D'Agostini's method (D'Agostini 2005J and the dashed line is the best 
fit obtained by the weighted x 2 method. Marginalizing over b, the likelihood values for the parameters a, a are a — 1.57 and 
a = 0.32 respectively. As in the previous cases we recalibrate also the E Pt i - E\ so correlation using the SNela data. Once 
such correlations have been calibrated, we can now use them to compute the estimated and calibrated GRBs Hubble diagram, 
according to the relations: 




E iso (l + z) 

4:TVSb0l0 



(22) 
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Figure 9. Best fit curves, relatively to the fit of the Amati correlation, superimposed on the data: the solid and dashed lines refer to 
the results obtained with the Bayesian and Levemberg - Marquardt estimators respectively. 



fl(z) = 25 + 51og d L {z). 

It turns out that there are seven GRBs deviating from the fiducial distance modulus fj, more than 2a; moreover also for the 
Amati correlation the estimated and calibrated HDs are consistent, as shown in Figs.QI)] and 1111 Moreover, in order to test 
the reliability of the Amati correlation with respect to the other correlations, we compare the distance modulus fi(z) for the 
calibrated Hubble diagram constructed of GRBs dataset used for fitting the five 2D correlations illustrated above and the 
calibrated Hubble diagram constructed of GRBs dataset used for fitting the Amati correlation. However there are only 21 
GRBs that appear in both the Schaefer (Schaefer 2007) and Amati ( Amati et al. 2011!)) samples. It turns out that there is 
full consistency between these datasets (see Fig. I12p . and the resulting distances are strongly correlated with the Spearmans 
correlation p = 0.92. We observe that our sample includes 4 GRBs with z > 5, namely 060522 at z — 5.11, 050904 at z = 6.29, 
060116 at z — 6.60 and 090423 at z — 8.23. As is apparent from Fig. llOl there is a clear gap between the z > 6 GRBs and the 
rest of the sample so that one could wonder whether these objects share the same properties with the other GRBs. Actually, 
even if extreme in redshift, all of them follow quite well the considered here correlation. In particular the observed properties 
of 090423 are in a good agreement with the Amati correlation. Therefore, we do not expect that they introduce a bias in the 
fiducial HD. Moreover, since they are at such a high redshift, they are not used in the calibration based on the local regression 
analysis. Nevertheless, their distance moduli in the fiducial and calibrated HDs are consistent. We therefore conclude that, 
notwithstanding the doubts on the validity of this 2D correlation at very high z, the inclusion of these three GRBs does not 
bias the results. At the end we use the calibrated GRBs Hubble diagram to test if our fiducial cosmological model is able to 
describe the background expansion up to redshifts z ~ 8. In the Bayesian approach to model testing, we explore the parameter 
space through the likelihood function : 

Cgrb (p) oc exp [~Xgrb (p)/2] , (23) 
with 



J^GRB 

Xgrb(p) = 

i=l 



GRB 



(24) 



where ogrb takes into account the intrinsic scatter inherited from the scatter of GRBs around the E Pi i - E{ ao correlation, and 

t-0.07 
-0.05' 



p denotes the set of model parameters - Ho, h in our case. The likelihood values are H = 0.95 ± 0.03 h = 0.62ij5;2I, which 



correspond to Q m — 0.26^q'q4 an d io(0) = —0.77^^. If we marginalize over h we obtain Ho — 0.93+ O 'o|, which corresponds 
to n m = 0.3ll°;°* and w(0) = -0.78+°'°*. 



4 DISCUSSION AND CONCLUSIONS 

Recently several interesting correlations among the Gamma Ray Burst (GRB) observables have been identified. Proper 
evaluation and calibration of these correlations are needed to use the GRBs as standard candles constraining the expansion 
history of the universe up to redshifts of z ~ 8. Here we used the GRB data set recently compiled by Schaefer (2007) to 
investigate, in a quintessential cosmological scenario, the fitting of six 2D correlations: namely Ti ag (time lag), V (variability), 
Speak (peak of the uF v spectrum), and trt (minimum rise time). In particular we have investigated the dependence of the 
calibration, on the assumed cosmological model and a possible evolution of the correlation parameters with the redshift. To 
check the robustness of the discovered correlations between some of the observed parameters of GRBs and their basic intrinsic 
properties, we have first tested the stability of values of the correlation parameters by considering an ad hoc definition of 
luminosity distance that gives much larger distances to objects at z > 2 than either of the fiducial models. It is apparent that 
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Figure 10. The calibrated GRBs Hubble diagram with overplotted the distance modulus predicted by the fiducial model. The full boxes 
correspond to the data set of 70 GRBs compiled by Amati et al. 2008, while the empty ones correspond to the added sample of 20 GRBs 
compiled by Amati et al. 2009. 



MSN 

Figure 11. Comparison of the distance modulus fi(z) for the calibrated and estimated GRBs Hubble diagram made up fitting the Amati 
correlation. 

our artificial luminosity distance is changing the values of the correlation parameters, however the difference is not dramatic 
at lcr of confidence, but becomes significant at higher levels. We then investigated the effects of cosmological parameters on 
the calibration parameters (a, b, <7i„t) when we explore the parameters-space of our fiducial cosmological model. It turns out 
that they very weakly depend on the assumed values of the cosmological parameters. Moreover our analysis shows only a 
modest preference for the E 1 -E p , L-Ti ag and Lx-T a correlations. Taken at face values, the maximum likelihood estimates 
for Oi n t favour the - E p correlation, which seems the most robust, but when one takes into account the 68% and 95% 
confidence ranges all the correlations overlap quite well. In our analysis we have assumed that the correlation parameters do 
not change with the redshift, which indeed spans a quite large range (from z — 0.125 up to z — 6.6). The limited number 
of GRBs prevents detailed exploration of the validity of this usually adopted working hypothesis, which we tested somewhat 
by investigating if the residuals correlate with the reshift. Moreover we tested the correlation E-,-E p with respect to the 
evolution with redshift, separating the GRB samples into four groups corresponding to the following redshift bins: z 6 [0, 1], 
z £ [1,2], z 6 [2,3] and z £ [3,7]. We thus maximized the likelihood in each bin of redshift and determined the best fit 
calibration parameters a, b with la errors and the intrinsic dispersion (Tint- We have not found any statistical evidence of a 
dependence of the (a,b,ai„t) parameters on the redshift, as recently also reported by Cardone et al. (2008) and Basilakos 
& Perivolaropoulos (2008). Once the six correlations have been fitted, we constructed the GRBs Hubble diagram, taking a 
weighted average of the distance modulus provided by each of the six 2D correlations considered. Although the above analysis 
has shown that the choice of the underlying cosmological model has only a modest impact on the final estimate of the distance 
modulus, we apply a local regression technique to estimate, in a model independent way, the distance modulus from the 
most updated SNIa sample, in order to calibrate the GRBs 2D correlations, including only GRBs with z < 1.4. The derived 
calibration parameters have been used to construct a new GRBs Hubble diagram, which we call the calibrated GRBs HD. We 
finally compare the estimated and calibrated GRBs HDs. It turns out that for the common GRBs they are fully statistically 
consistent thus indicating that both of them are not affected by any systematic bias induced by the different fitting procedures. 
This means that the high redshift GRBs can be used to test different models of dark energy. We also constructed the Hubble 
diagram using the E Pt i - E[ so correlation, which, recently has been shown to be useful in cosmological considerations. In our 
analysis we used the sample of 95 long GRB/XRF only, compiled in (Amati et al. 2008) and (Amati et al. 2009), which covers 
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Figure 12. Comparison of the distance modulus /J,(z) for the calibrated GRBs Hubble diagram made up by fitting the five 2 — D 
correlations illustrated above, nsNold, and the estimated GRBs Hubble diagram made up by fitting the Amati correlation, fii ao , limited 
to the 21 GRBs that appear in both the Schaefer and Amati samples (Schaefer 2007, Amati et al. 2008, Amati et al. 2009). It turns out 
that these datasets are fully consistent. 

a broad range of z, and includes GRB 090423, the new high-z record holder of Gamma-ray bursts. As above we calibrated 
this correlation maximizing our likelihood function. It turns out that our likelihood method gives a — 1.56, b = 52.77 and 
Oint = 0.31. Applying the same Bayesian technique we recalibrate also the E Pi i - Ei so correlation using the SNIa data. It 
turns out that there are seven GRBs deviating from the fiducial distance modulus fi more than 2a; so that also for the Amati 
correlation the fiducial and calibrated HDs are consistent. We finally use the calibrated GRBs Hubble diagram to test if our 
fiducial quintessential cosmological model is able to describe the background expansion up to redshifts z ~ 8. The likelihood 
values correspond to Q m — 0.261q'q4 an d w(0) = — 0.74q"q3 01 . In conclusion one can say that the correlation parameters that 
relate some of the observed parameters of GRBs with their basic intrinsic parameters only weakly depend on the background 
cosmological model. This in turn means that the Gamma Ray Bursts could be used to derive values of the basic cosmological 
parameters and in the future with much larger data set it should be possible to get more information about the nature of dark 
energy. Since the discovered correlations only weakly depend on the background cosmological model they should be created 
during the burst itself and therefore they also provide new constrains on models of the Gamma Ray Burst. 
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Table 1. Calibration parameters (a, b), intrinsic scatter <Tj n t, median, 5 me d, and root mean square of the best fit residuals, S rrns for the 
2D correlations logR = a log Q + b evaluated in two fiducial cosmological models (the ACDM and the quintessential exponential scalar 
field models) and in our crazy model. Columns are as follows : 1. id of the correlation; 2. number of GRBs used; 3. maximum likelihood 
parameters; 4,5 median value and 68% and 95% confidence ranges for the parameters (a, <Ji„t); 6,7 median and root mean square of the 
residuals. 
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